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Abstract 

We report the first study of restoration of rotational symmetry and fluctuations of the quadrupole 
deformation in the framework of relativistic mean-field models. A model is developed which uses 
the generator coordinate method to perform configuration mixing calculations of angular momen- 
tum projected wave functions, calculated in a relativistic point-coupling model. The geometry is 
restricted to axially symmetric shapes, and the intrinsic wave functions are generated from the 
solutions of the constrained relativistic mean-field + BCS equations in an axially deformed oscil- 
lator basis. A number of illustrative calculations are performed for the nuclei 194 Hg and 32 Mg, in 
comparison with results obtained in non-relativistic models based on Skyrme and Gogny effective 
interactions. 

PACS numbers: 21.60.Jz, 21.10.Pc, 21.10.Re, 21.30.Fe 
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I. INTRODUCTION 



The rich variety of nuclear shapes far from stability has been the subject of extensive 
experimental and theoretical studies. The variation of ground-state shapes in an isotopic 
chain, for instance, is governed by the evolution of shell structure. In particular, far from the 
/^-stability line the energy spacings between single-particle levels change considerably with 
the number of neutrons and/or protons. This can result in reduced spherical shell gaps, 
modifications of shell structure, and in some cases spherical magic numbers may disappear. 
For example, in neutron-rich nuclei N = 6,16,34... can become magic numbers, whereas 
iV = 8,20,28... disappear. The reduction of a spherical shell closure is associated with the 
occurrence of deformed ground states and, in a number of cases, with the phenomenon of 
shape coexistence. 

Both the global shell-model approach and self-consistent mean-field models have been 
employed in the description of shell evolution far from stability. The basic advantages of the 
shell model is the ability to describe simultaneously all spectroscopic properties of low-lying 
states for a large domain of nuclei, effective interactions that can be related to two- and 
three-nucleon bare forces, and a description of collective properties in the laboratory frame. 
On the other hand, since effective interactions strongly depend on the choice of active shells 
and truncation schemes, there is no universal shell-model interaction that can be used for all 
nuclei. Moreover, because a large number of two-body matrix elements has to be adjusted 
to data, extrapolations to exotic systems far from stability cannot be very reliable. Heavy 
exotic nuclei with very large valence spaces require calculations with matrix dimensions that 
are far beyond the limits of current shell model variants. 

Properties of heavy nuclei with a large number of active valence nucleons are best de- 
scribed in the framework of self-consistent mean-field models. A variety of structure phe- 
nomena, not only in medium-heavy and heavy stable nuclei, but also in regions of exotic 
nuclei far from the line of /3-stability and close to the nucleon drip-lines, have been suc- 
cessfully described with mean-field models based on the Gogny interaction, the Skyrme 
energy functional, and the relativistic meson-exchange effective Lagrangian jjIQ- The self- 
consistent mean-field approach to nuclear structure represents an approximate implemen- 
tation of Kohn-Sham density functional theory, which enables a description of the nuclear 
many-body problem in terms of a universal energy density functional. This framework, 
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extended to take into account the most important correlations, provides a detailed micro- 
scopic description of structure phenomena associated with shell evolution in exotic nuclei. 
When compared to the shell model, important advantages of the mean-field approach include 
the use of global effective nuclear interactions, the treatment of arbitrarily heavy systems 
including superheavy elements, and the intuitive picture of intrinsic shapes. 

A quantitative description of shell evolution, in particular the treatment of shape coex- 
istence phenomena, necessitates the inclusion of many-body correlations beyond the mean- 
field approximation. The starting point is usually a constrained Hartree-Fock plus BCS 
(HFBCS), or Hartree-Fock-Bogoliubov (HFB) calculation of the potential energy surface 
with the mass quadrupole components as constrained quantities. In most applications cal- 
culations have been restricted to axially symmetric, parity conserving configurations. The 
erosion of spherical shell-closures in nuclei far from stability leads to deformed intrinsic states 
and, in some cases, mean-field potential energy surfaces with almost degenerate prolate and 
oblate minima. To describe nuclei with soft potential energy surfaces and/or small energy 
differences between coexisting minima, it is necessary to explicitly consider correlation ef- 
fects beyond the mean-field level. The rotational energy correction, i.e. the energy gained by 
the restoration of rotational symmetry, is proportional to the quadrupole deformation of the 
intrinsic state and can reach several MeV for a well deformed configuration. Fluctuations of 
quadrupole deformation also contribute to the correlation energy. Both types of correlations 
can be included simultaneously by mixing angular momentum projected states correspond- 
ing to different quadrupole moments. The most effective approach for configuration mixing 
calculations is the generator coordinate method (GCM), with multipole moments used as 
coordinates that generate the intrinsic wave functions. 

In a series of rece. papers gflflQ , the angular momentum projected GCM with 
the axial quadrupole moment as the generating coordinate, and intrinsic configurations 
calculated in the HFB model with the finite range Gogny interaction, has been applied in 
studies of shape-coexistence phenomena that result from the erosion of the N = 20 and 
N = 28 spherical shells in neutron-rich nuclei. Good agreement with experimental data 
has been obtained for the 2 + excitation energies, and B(E2) transition probabilities of the 
N = 28 neutron- rich isotones . The systematic study of the ground and low- lying excited 
states of the even-even 20_40 Mg 4] is particularly interesting, because this chain of isotopes 
includes three spherical magic numbers N = 8, 20, 28. It has been shown that the N = 8 shell 
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closure is preserved, whereas deformed ground states are calculated for N = 20 and N = 28. 
In particular, the ground state of 32 Mg becomes deformed as a result of a fine balance 
between the energy correction associated with the restoration of rotational symmetry and 
the correlations induced by quadrupole fluctuations. In a similar analysis of the chain of 
even-even isotopes 20-34 Ne p], it has been shown that the ground state of the N = 20 
nucleus 30 Ne is deformed, but less than the ground state of its isotone 32 Mg. The model 
has recently been applied in an analysis of shape coexistence and quadrupole collectivity in 
the neutron- deficient Pb isotopes A good qualitative agreement with available data has 
been found, especially for rotational bands built on coexisting low-lying oblate and prolate 
states. 

r-i 

Another very sophisticated model [7| which extends the self-consistent mean-field ap- 
proach by including correlations, is based on constrained HF+BCS calculations with Skyrme 
effective interactions in the particle-hole channel and a density-dependent contact force in 
the pairing channel. Particle number and rotational symmetry are restored by projecting 
self-consistent mean-field wave functions on the correct numbers of neutrons and protons, 
and on angular momentum. Finally, a mixing of the projected wave functions corresponding 
to different quadrupole moments is performed with a discretized version of the generator co- 
ordinate method. The model has been successfully tested in the study of shape coexistence 
in 16 0, and in the analysis of the coexistence of spherical, deformed, and superdeformed 
states in 32 S, 36 Ar, 38 Ar and 40 Ca jj. For the doubly-magic nucleus 16 this parameter-free 
approach provides a very good description of those low-spin states which correspond to axi- 
ally and reflection-symmetric shapes, and allows the interpretation of their structure in terms 



of self-consistent np — nh states. A very important recent application is the study of 
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lying collective excitation spectra of the neutron-deficient lead isotopes 182_194 Pb 
A configuration mixing of angular-momentum and particle-number projected self-consistent 
mean-field states, calculated with the Skyrme SLy6 effective interaction, qualitatively re- 
produces the coexistence of spherical, oblate, prolate, and superdeformed prolate structures 
in neutron-deficient Pb nuclei. 

Even though the self-consistent relativistic mean-field (RMF) framework has been em- 
ployed in many studies of deformed nuclei, applications of meson-exchange and point- 
coupling models have so far been restricted to the mean-field level. In this work we report 
the first study of restoration of rotational symmetry and fluctuations of the quadrupole 
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deformation in the framework of relativistic mean-field models. We perform a GCM con- 
figuration mixing of angular momentum projected wave functions that are calculated in a 
relativistic point-coupling model. 

In Section |n] we present an outline of the relativistic point-coupling model which will be 
used to generate mean-field wave functions with axial symmetry, introduce the formalism 
of the generator coordinate method, and describe in detail the procedure of configuration 
mixing of angular momentum projected wave functions. In Section ITTT1 our model for GCM 
configuration mixing is investigated in a study of quadrupole dynamics in the nucleus 194 Hg, 
and 32 Mg is used as a test case for the configuration mixing calculation of angular momentum 
projected states. Section |TV| summarizes the results of the present investigation and ends 
with an outlook for future studies. 



II. CONFIGURATION MIXING OF ANGULAR MOMENTUM PROJECTED 
MEAN-FIELD WAVE FUNCTIONS 

In this section we review the self-consistent relativistic point-coupling model which will 
be used to generate constrained mean-field states, and the solution of the corresponding 
single-nucleon Dirac equation in an axially symmetric harmonic oscillator basis. Starting 
with a short outline of the generator coordinate method, we describe the technical details 
of the configuration mixing of angular momentum projected wave functions. 



A. The relativistic point-coupling model 

Most applications of the self-consistent relativistic mean-field framework have used the 
finite-range meson-exchange representation, in which the nucleus is described as a system of 
Dirac nucleons coupled to exchange mesons and the electromagnetic field through an effective 
Lagrangian. A medium dependence of the effective nuclear interaction can be introduced ei- 
ther by including non-linear meson self-interaction terms in the Lagrangian, or by assuming 
an explicit density dependence for the meson- nucleon couplings Q|. An alternative represen- 
ation is formulated in terms of point-coupling (PC) (contact) nucleon-nucleon interactions 
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161 ] . In RMF-PC models the medium dependence of the interaction can be 



taken into account by the inclusion of higher order interaction terms, for instance six-nucleon 
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vertices (?pip) 3 , and eight-nucleon vertices (tptp) 4 and [( , ijj^ lx ip)('ip'j^ip)} 2 , or it can be encoded 
in the effective couplings, i.e. in the strength parameters of the interaction in the isoscalar 
and isovector channels. When employed in studies of ground-state properties of finite nuclei, 
the two representations produce results of comparable quality. The point-coupling formu- 
lation, however, avoids some of the constraints imposed in the meson-exchange picture as, 
for instance, the use of the fictitious sigma-meson in the isoscalar-scalar channel. The self- 
consistent PC models are also closer in spirit to the nuclear density functional theory, in 
which the exact energy functional, including higher-order correlations, is approximated with 
powers and gradients of ground-state nucleon densities. The point-coupling representation, 
with medium-dependent vertex functions, provides a natural framework in which chiral ef- 
fective field theory can be employed to construct the nuclear energy density functional, thus 
establishing a link between the rich nuclear phenomenology and the underlying microscopic 
theory of low-energy QCD [l^Ql- 

A detailed description of the point- coupling model that we use in this work can be found, 
for instance, in Ref. together with a thorough discussion of the choice of various 

parameter sets that determine the effective interactions. Here we only outline the essential 
features of the model and of its mean-field solution for a deformed axially symmetric nucleus. 

The relativistic point- coupling Lagrangian is built from basic densities and currents bi- 
linear in the Dirac spinor field ip of the nucleon: 

tpOrTi/j , T e{l,Ti} , r G {1,7,*, 75, 757/*, <v> • (!) 

Here Tj are the isospin Pauli matrices and T generically denotes the Dirac matrices. The 
interaction terms of the Lagrangian are products of these bilinears. Although a general 
effective Lagrangian can be written as a power series in the currents ipOrTip and their 
derivatives, it is well known from numerous applications of relativistic mean-field models 
that properties of symmetric and asymmetric nuclear matter, as well as empirical ground 
state properties of finite nuclei, constrain only the isoscalar-scalar (S), the isoscalar- vector 
(V), the isovector- vector (TV), and to a certain extent the isovector-scalar (TS) channels. 
In this work we consider a model with four-, six-, and eight-fermion point couplings (contact 
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[16], 



(2) 



interactions) |l6j], denned by the Lagrangian density: 
£ = C bee + £ 4f + £ hot + £ der + £ em , 

£ free = ^(yy^-m)^, 

£hot = _| /5s( ^ ) 3_ | TS (^) 4 - J7V[(^)(^)] 2 , 

£ em = -eA^[(l - r 3 )/2] 7 ^ - ± F^F^. 

Vectors in isospin space are denoted by arrows, and bold-faced symbols will indicate vectors 
in ordinary three-dimensional space. In addition to the free nucleon Lagrangian jCf ree , the 
four-fermion interaction terms contained in £ 4 f , and higher order terms in £hot> when applied 
to finite nuclei the model must include the coupling £ em of the protons to the electromagnetic 
field A^, and derivative terms contained in £der- In the terms d^^Tip) the derivative is 
understood to act on both %p and if). One could, of course, construct many more higher 
order interaction terms, or derivative terms of higher order, but in practice only a relatively 
small set of free parameters can be adjusted from the data set of ground state nuclear 
properties. 

The single-nucleon Dirac equation is derived from the variation of the Lagrangian (J2J) 
with respect to ip 

{ a [_ zV - V(r)] + V(r) + (3(m + S(r))} ^(r) = e^r) . (3) 
The scalar and vector potentials 

S(r) = Z s (r) + T 3 Z TS (r) , (4) 

V^r) = ^(r) + r 3 ^ TV (r), (5) 

contain the nucleon isoscalar-scalar, isovector-scalar, isoscalar-vector and isovector-vector 
self-energies defined by the following relations: 

£s = asps + PsPs + Ispl + $s&ps , (6) 



T, TS = cutsPts + Sts^Pts , (7) 
W = a v f + -yvUvfW + SvAj* ~ e ^^^ > (8) 
££y = a Tv3 % + 5tv^3tv > ( 9 ) 

respectively. Because of charge conservation, only the 3 — rd component of the isovector 
densities and currents contributes to the nucleon self-energies. The local densities and 
currents are calculated in the no-sea approximation 

A 

Ps(r) = Y,4(r)mr) , (10) 
i=i 

A 

Pts(t) = X>J(r)/3r 3 ^(r) , (11) 
i=i 

^(r) = 5>J(r)/Wi(r) , (12) 

i=l 
A 

Jrvir) = 52l>kr)Pfnfil>i(r) ■ (13) 
i=i 

For a nucleus with A nucleons, the summation runs over all occupied states in the Fermi 
sea, i.e. only occupied single-nucleon states with positive energy explicitly contribute to 
the nucleon self-energies. The energy momentum tensor determines the total energy of the 
nuclear system 

Ermf = J dr £ RM F{r) 

= ^ / dr ipj{r) (— ictV + f3m) ipi(r) 
i J 

+ /* (f A + f pi + f ^ + 5 i" sA " s + + f'w"' 2 + s i^ Af 

+ —JtvUtv),* + —j^ v A(j TV )^ + —Pts + —Pts^Pts + -p P A J , (14) 

where p p denotes the proton density, and A is the Coulomb potential. 

In this work we only consider even-even nuclei that can be described by axially symmetric 
shapes. It is therefore convenient to work in cylindrical coordinates 

x = r^cos0, y = r^sin0 and z. (15) 

In addition to axial symmetry, parity, symmetry with respect to the operator e~ lnJy , and 
time-reversal invariance are imposed as self-consistent symmetries. Time-reversal invariance 
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implies that spatial components of the currents vanish in the nuclear ground state. The 
resulting single-nucleon Dirac equation reads 



{-ictV + V(r) + (3(m + S{r))}^{r) = e^r) . 



(16) 



The eigensolutions are characterized by the projection of the total angular momentum along 
the symmetry axis (£1%), the parity (tt^), and the z-component of the isospin (tj). The Dirac 
spinor has the following form 

I /+(^,r ± )e i ^- 1 /2)* ^ 
/-(^,r ± )e^ +1 / 2 ^ 

^(^-1/2)0 



fi(r,s,t) 



m(r,s,t) 
For each solution with positive Q 



'2tx 



i(Qi+l/2)(t> 



xtM 



(17) 



the corresponding degenerate time-reversed state 



■01 — Ttpi — {-fr, /+, gr, - g +- -a} 



(18) 



(19) 



is obtained by acting with the time-reversal operator T = ia y K. For even-even nuclei, the 
time-reversed states i and % have identical occupation probabilities. 

The single-nucleon Dirac eigenvalue equation is solved by expanding the spinors fa and_ 
Eq. (|T7|) in terms of eigenf unctions of an axially symmetric harmonic oscillator potential 



V osc (z,r ± ) = l -Mu? z z 2 + X -Muj\r\ . 



(20) 



Imposing volume conservation, the two oscillator frequencies hu z and hu± can be expressed 
in terms of the deformation parameter f3 and the oscillator frequency ItjjJq 



and 



(21) 



The corresponding oscillator length parameters are b z = a/ h/Mu z and b± = \J h/Mu±. 
Because of volume conservation, b\b z = 6jj. The basis is now specified by the two constants 
hjjQ and flo, and basis states are characterized by the set of quantum numbers 



| a >= \n z ,n±,A,m s > , 



(22) 



where n z and n± denote the number of nodes in the z and r± directions, respectively. A and 
m s are the components of the orbital angular momentum and the spin along the symmetry 
axis. The eigenvalue of j z - the ^-projection of the total single-nucleon angular momentum 

tt = A + m s , (23) 

and the parity is given by 

7T = (-1)»*+ A . (24) 
The eigenfunctions of the axially symmetric harmonic oscillator potential read 

AT /V A 1 

*«(r, s) = %^(0^ 2/2 ^^ A / 2 L A x (^) e ^/ 2 -Le^ Xms ( S ) , (25) 



b ± 1 n±yu 2vr 
with ^ = z/b z and rj = r\jb\. The Hermite polynomials H n (£), and the associated Laguerre 
polynomials L^(rj), are defined in Ref. j^J. The normalization factors are given by 



N n = ; 1 and = J , n±l - . (26) 



The large and small components of the single-nucleon Dirac spinor Eq. ()17|) are expanded 
in terms of the eigenfunctions Eq. (|2*5T) 

1 / f + (z r i ) e i ^~ 1 / 2 ^ \ amax 

fi{r,*,t) = -= r; ' . (n+l/2U = E /r$ a (r,,) Xtl (t) , (27) 

In order to avoid the onset of spurious states, the quantum numbers a max and a max are 
chosen in such a way that the corresponding major oscillator quantum numbers N = n z + 
2n± + A are not larger than N s h + 1 for the expansion of the small components, and not 
larger than N s h for the expansion of the large components jioj ]. 

For an axially deformed nucleus the map of the energy surface as a function of the 
quadrupole moment is obtained by imposing a constraint on the mass quadrupole moment. 
The method of quadratic constraint uses an unrestricted variation of the function 

<H> + ^(<Q>- q y , (29) 
where < H > is the total energy, < Q > denotes the expectation value of the mass 



quadrupole operator, q is the deformation parameter, and C is the stiffness constant 
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21] 



In addition to the self-consistent mean-field potential, for open-shell nuclei pairing cor- 
relations have to be included in the energy functional. In this work we do not consider 
nuclear systems very far from the valley of /^-stability, and therefore a good approximation 
for the treatment of pairing correlations is provided by the BCS formalism. Following the 
prescription from Ref. [16j] . we use a 6 force in the pairing channel, supplemented with a 
smooth cut-off determined by a Fermi function in the single-particle energies. The pairing 
contribution to the total energy is given by 

£& = j &r)dr = ^ff K ; in) (r)K p(n) (r)dr , (30) 

for protons and neutrons, respectively. K p ^(r) denotes the local part of the pairing tensor, 
and V p (n) is the pairing strength parameter. Of course, for open-shell nuclei the expressions 
Eqs. (110)1 - (fT3*|) f° r the local densities and currents include the occupation factors of single- 
nucleon states. Finally, the expression for the total energy reads 

E tot = f [£ RMF (r) + £ p patr (r) + £; air (r)] dr . (31) 

The center-of-mass correction has been included by adding the expectation value 

(P 2 ) 

to the total energy. P cm is the total momentum of a nucleus with A nucleons. 



B. The generator coordinate method 

The generator coordinate method (GCM) is based on the assumption that, starting from 
a set of mean-field states \<j>(q)) which depend on a collective coordinate q, one can build 
approximate eigenstates of the nuclear Hamiltonian 

l*«> = X>(<&) M&)> ■ (33) 
j 

A detailed review of the GCM can be found in Chapter 10 of Ref. [2l|. In this work the 
basis states \4>(q)) are Slater determinants of single- nucleon states generated by solving the 
constrained relativistic mean-field + BCS equations, as described in the previous section. 
This means that we use the mass quadrupole moment as the generating coordinate q. The 
axially deformed mean-field breaks rotational symmetry, so that the basis states \4>{q)) are 
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not eigenstates of the total angular momentum. Of course, in order to be able to compare 
theoretical predictions with data, it is necessary to construct states with good angular 
momentum 

where P(j K denotes the angular momentum projection operator 

P mk = J dQD{* K (Q)R(Q) . (35) 

Integration is performed over the three Euler angles a, ft, and 7. D J MK {Vt) = 
e~ iMa d J MK {ft)e~ iK ' 1 is the Wigner function 22], and R(Q) = e~ iaSz e'^ e~^ jz is the rotation 
operator. The weight functions fi K {qj) are determined from a variational calculation, 





H | 




(H M \ 


a 1 



SE J = i M',';/ = , (36) 



i.e. by requiring that the expectation value of the energy is stationary with respect to an 
arbitrary variation Sf£ K ■ This leads to the Hill- Wheeler equation 

E fi K (%) {^)\HP J mk I0(?i)> - K (<P{ qi )\P J MK I0(?i)>) = • (37) 

The restriction to axially symmetric configurations (J z \(f>(q)) = 0) simplifies the problem 
considerably, because in this case the integrals over the Euler angles a and 7 can be per- 
formed analytically. For an arbitrary multipole operator Qx^ one thus finds 

2 J + 1 f n - 

(<P(q t )\ Q^PliK I0fe)> = — g— Sm-^kq J sin/3^; (/3) (0(%)| Qa^"^ \<j>{ qj )) dft . 

(38) 

By using the identity e l/3Jy = e~ lnJz e~' ll3Jy e' lnJz , together with parity, and the symmetry with 
respect to the operator e~ mJy , the integration interval in Eq. (|3*Hj) can be reduced from [0, n] 
to [0,7r/2] 

{^Qx^mM^)) = (2J+1) 1 + { ~ 1)J 5 M ^S K0 

^ smftdt^ift) (^IQx^y dft . (39) 

We notice that this expression vanishes for odd values of angular momentum J, i.e., the 
projected quantities are defined only for even values of J. 



12 



The norm overlap kernel 



Af J (q t ,q,) = (<P{qi)\P J MK \<P{q 3 )) 

1 + 



(27 i i ' ' v o l)J 5 M0 6 K0 [ ' sm(3d J *{[3) (<P{ qi )\e-^ \<P{ qj )) d{3 , (40) 
1 Jo 



can be evaluated by employing the generalized Wick theorem 



E, 23, 24, 3 



n(q i: qf, (3) = (4>{qi)\e \<j)( gj )) = ±Jdet M ab (q u q 3 \ 0) 



(41) 



The overlap matrix is defined: 



Nab(qi, qj] P) = u a (qi)R ab (qi, qj] (3)u b (q 3 ) + v a (qi)R ab (q u q 5 ] /3)v b (qj) , 
where u and v denote the BCS occupation probabilities, and the matrix R reads 



(42) 



Rab{qi,qj;P) = J ^ a {r;q l )e i/3J ^(r; qj)dr . (43) 

If the expansions Eqs. (J27|) and (J28j) are inserted in the expression above, the evaluation of 
the matrix R reduces to the calculation of the matrix elements of the rotation operator in 
the basis of the axially symmetric harmonic oscillator 



Rab(qi,qj) = J2fa(qi)f b (qj 

a,/3 



a 



(a 



-i(3J y 



(44) 



The simplest way to evaluate these matrix elements is to express the eigenf unctions of the 
axially symmetric harmonic oscillator in the spherically symmetric oscillator basis. The 
transformation from the spherical to the axially deformed basis is given by the following 
expression 



|fi Q A a n"n°) 



1 Q r 

/ j °0. a K a n a 1 n\ 
nlj 



\nljQ a ) , 



(45) 



with the transformation coefficients S^ AaTiana given in Ref. 



It must be emphasized 

that this transformation is only possible if uo z = u± (i.e. (3q = in Eq. (|21|) ) [^Sl- I n 
addition, the same oscillator frequency hu> has to be used for each value of the generating 
coordinate q in order to avoid completeness problems in the GCM calculations j^. We 
have used Jtluq = 41A -1 / 3 . Since the choice of the two basis parameters ftiuo and (3q cannot 
be optimized, the convergence of the results should be carefully checked as a function of the 
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number of oscillator shells used in the expansion of the Dirac spinors. The expression for 
the matrix elements is simply 



(a e 



-i/3J y 



Eqnlj qnlj ij ,a\ 

nlj 



(46) 



where d 3 ^ aVL ((3) denotes the Wigner rotation matrix (22 |. Adore general transformation coef- 



ficients for the case uj± ^ u z have been derived in Ref. 
and have not been used in the present analysis. 
The Hamiltonian kernel 



29j, but they are rather complicated 



(2J + l) 1 + ( 2 1)J 5 M0 5 K o sin (3d J *(f3) (0(g,)| He~^ |#?;)) d/3 , (47) 



can be calculated from the mean-field energy functional Eq. (JTlJl 
the modified densities 



23 



24123, 

provided 



Ps(r;qi,qj,/3 
pTs(r;q i} q j} f3 

jTv(r;qi,qj,/3 



v Mi)vbi.qj W ba 1 (q i} qf, (3)^l(r; g<)(-iaV + mf3)e i/3J ^"ip b (r; q 5 

a.b 

M^MVjWba 1 ^ Qj\ P)fl>l(r; q^pe-^Mr; qj) , (49) 

a.b 

v a{qi)vb{qj)M^{qi, qf, (3)^ a (r; qi )f3T^ J ^ b {r- q,) , (50) 

a.b 

VaiQiMqjW^iqi, Qj ; /3)^(r; q^^Hb^ Qj ) , (51) 

a.b 

VaiqiHiqjWiaiQu PWi(r; q l )Pl fl r 3 e- i ^Mr; q 3 ) , (52) 



a.b 



are used when evaluating the expression 

h(q h q 3 ;P) = (<f>( qi )\ H \<f>( qj )) = J £ tot (r; q h Qj , (3)dr . (53) 

The computational task of evaluating the Hamiltonian and norm overlap kernels can be 
reduced significantly if one realizes that states with very small occupation probabilities give 
negligible contributions to the kernels. Such states can be excluded from the calculation, 
a Ildtll e d e tai ,of tlli8pr o C e dur e Canb efou„d in Re f ,SQ. 

An additional problem arises from the fact that the basis states \4>{qj)) are not eigenstates 
of the proton and neutron number operators Z and N. The adjustment of the Fermi energies 
in a BCS calculation ensures only that the average value of the nucleon number operators 
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corresponds to the actual number of nucleons. Consequently, the wave functions |^/^ M ) 
are generally not eigenstates of the nucleon number operators and, moreover, the average 
values of the nucleon number operators are not necessarily equal to the number of nucleons 
in a given nucleus. This happens because the binding energy increases with the average 
number of nucleons and, therfore, an unconstrained variation of the weight functions in 
a GCM calculation will generate a ground state with the average number of protons and 
neutrons larger than the actual values in a given nucleus. In order to restore the correct 



mean values of the nucleon numbers, we follow the usual prescription 
the Hill- Wheeler equation by replacing h(qi, qf (3) with 
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and modify 



h'(qi, qj; (3) = h(qi, qf (3) - X p [z(q i} qj; (3) - z ) - A n [n(q h q,-; (3) - n ] , (54) 

where 

z(q l ,q j] p) = ( ( f ) (q l )\Ze-^ } y\4 > (q,)) and n(fc Qj ; (3) = (0(ft)| Ne^ |0(<&)> . (55) 

A p ( n ) is the proton (neutron) Fermi energy, while Zq and n denote the desired number of 
protons and neutrons, respectively. 
The Hill- Wheeler equation 

3 3 

presents a generalized eigenvalue problem, and thus the weight functions fi{q%) are not 
orthogonal and cannot be interpreted as collective wave functions for the variable q. It is 
useful to re-express Eq. (jHHj) in terms of another set of functions, g^iQi), defined by 

3 

With this transformation the Hill- Wheeler equation defines an ordinary eigenvalue problem 

52n J (qi,<a)9t(qj) = E a g J M^ , (58) 

3 

with 

n J (q r , qj ) = ^(^)- 1 /2 te;gfc )^ (gfc)g/ ) ( ^)-V 2(ghgj) . (59) 

k,l 

The functions g^iQi) are orthonormal and play the role of collective wave functions. 
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In practice, the first step in the solution of Eq. (|56|) is the diagonalization of the norm 
overlap kernel M J {qi,qj) Eq. (jJDj) 

S ^N' J (qi,qj)u k (q j ) = n k u k (qi) . (60) 

3 

Since the basis functions \4>{qi)) are not linearly independent, many of the eigenvalues n k are 
very close to zero. The corresponding eigenf unctions u k (qi) are rapidly oscillating and carry 
very little physical information. However, due to numerical uncertainties, their contribution 
to Ti, J (qi,qj) can be large, and these states should be removed from the basis. From the 
remaining states one builds the collective Hamiltonian 

h£ = t- * I>fc(g^%,giH(gj) , (6i) 

which is subsequently diagonalized 

= Eigf . (62) 

k,l 

The solution determines both the ground state energy, and the energies of excited states, 
for each value of the angular momentum J. The collective wave functions g^iq), and the 
weight functions are calculated from the norm overlap eigenf unctions 

g J a ( qi ) = J29i a M<n), (63) 
l 

and 

Once the weight functions /„(<?) are known, it is straightforward to calculate all physical 
observables, such as transition probabilities and spectroscopic quadrupole moments |4(. The 
reduced transition probability for a transition between an initial state (J^otj), and a final 
state (Jf,aif), reads 

2 



B{E2;J i a i -> J /a/ ) 



e 2 



E /«/*(?/) IQ2I fe) /£(* 



2J, + 1 

and the spectroscopic quadrupole moment for a state (Jo) is defined 



(65) 



<T ec (J,«) = eJ 1 -^ l J2J \ W IQ2I |Jfc> (<&) • (66) 



J -J 
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Since these quantities are calculated in full configuration space, there is no need to introduce 
effective charges, hence e denotes the bare value of the proton charge. In order to evaluate 
transition probabilities and spectroscopic quadrupole moments, we will need the reduced 
matrix element of the quadrupole operator 



r\in(5d%M {<P(q s )\Q2,e-^y |0(g,)> . (67) 
Jo 

III. ILLUSTRATIVE CALCULATIONS 

In this section we perform several illustrative configuration mixing calculations that will 
test our implementation of the generator coordinate method, as well as the angular mo- 
mentum projection. The intrinsic wave functions that are used in the configuration mix- 
ing calculation have been obtained as solutions of the self-consistent relativistic mean-field 
equations, subject to constraint on the mass quadrupole moment. The interaction in the 
particle- hole channel is determined by the effective point-coupling Lagrangian Eq. (j2J), and a 
density-independent 5-force is used as the effective interaction in the particle-particle chan- 
nel. Pairing correlations are treated within the BCS framework. 

Among a number of self-consistent RMF-PC models that have been considered over the 
last ten years, a few reliable and accurate phenomenological parameterizations have been 
adjusted and applied in the description of ground state properties of finite nuclei on a quanti- 
tative level. In particular, based on an extensive multiparameter \ 2 minimization procedure, 
Biirvenich et al. have adjusted the PC-F1 set of couplin g co nstants for an effective point- 
coupling Lagrangian with higher order interaction terms [16]. While the Lagrangian of Eq. 
(J2J) contains 11 adjustable coupling constants, the PC-F1 effective interaction corresponds to 
a restricted set of 9 coupling parameters and does not include the isovector-scalar channel. 
In addition, the effective pairing interaction is determined by the strength parameters V p 
and V n , for protons and neutrons, respectively. The parameters in the particle-hole and 
particle-particle channels have been adjusted to ground state observables (binding ener gies , 
charge radii, diffraction radii, surface thickness, and pairing gaps) of 17 spherical nuclei [l6(. 

The PC-F1 interaction has been tested in the analysis of the equations of state of sym- 
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metric nuclear matter and neutron matter, binding energies and form-factors, and shell- 
structure-related ground-state properties of several isotopic and isotonic chains. This in- 
teraction has also been employed in relativistic quasiparticle random phase approximation 
calculations of multipole giant resonances 31]. A comparison with data has shown that the 
RMF-PC model with the PC-F1 interaction can compete with the best phenomenological 
finite-range meson-exchange interactions. It should be noted, however, that PC-F1 exhibits 
a relatively large volume asymmetry at saturation, resulting in a very stiff equation of state 
for neutron matter and too large values for the neutron skin in finite nuclei. Modern meson- 
exchange effective interactions, on the other hand, include an explicit medium dependence in 



both isoscalar and isovector channels 



3 



34j , and thus provide an improved description 



of asymmetric nuclear matter and neutron matter, and realistic values of the neutron skin 
in finite nuclei. 



A. Test of the generator coordinate method: 194 Hg 

Our first example is a test of the generator coordinate method in configuration mixing 
calculations for the nucleus 194 Hg. At this stage we do not consider angular momentum 
projection yet. The results of the test for the ground and excited states will be directly 
compared with the classical analysis of the GCM in the study of shape isomerism in 194 Hg 
by Bonche et al. j^J. It has to be emphasized, however, that the calculated GCM energies 
cannot be compared with data on a quantitative level, because without angular momentum 
projection not only the rotational energy correction is missing, but also the overlaps between 



states which belong to prolate and oblate minima are significantly reduced 35]. 

The GCM basis is constructed from self-consistent solution of the constrained single- 
nucleon Dirac equation on a regular mesh in the generating coordinate - the mass quadrupole 
moment: from q = —40b to q = 80b, with a spacing of Aq = 2b. The GCM basis thus 
consists of 61 intrinsic states. The large and small components of Dirac spinors are expanded 
in terms of the axially symmetric oscillator eigenfunctions. As already pointed out in Sec. 



Ill Bl the same oscillator frequency Hujq is used for each value of the generating coordinate 
q, and additionally the condition u z = uj± is imposed. Since the basis parameters are fixed 
to ftujQ = 41A -1 / 3 and (3q = 0, rather then optimized, the convergence of the results with 
respect to the number of major oscillator shells used in the expansions Eqs. ()27j) and (|2*K|) 
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has to be checked carefully. 

In Fig. ^ we display the binding energy curves for 194 Hg, as functions of the mass 
quadrupole moment, calculated by expanding the Dirac spinors in 10, 12, 14, 16 and 18 
oscillator shells. Obviously, at least 14 oscillator shells are necessary in order to obtain 
convergence for deformations smaller than q = 35b. Larger deformations require at least 16 
major oscillator shells. The absolute minimum of the binding energy curve corresponds to a 
slightly oblate shape (q = —10b). An additional shallow minimum at excitation energy «1.5 
MeV is found on the prolate side q = 6 b. At much larger deformation, q = 45 b, we find 
a third, superdeformed minimum 4.2 MeV above the first minimum of the binding energy 
curve. The deformations at which the three minima occur are in quantitative agreement 
with those calculated with the non-relativistic constrained Hartree-Fock plus BCS model of 
Ref. (2^, using the SIII Skyrme effective interaction. The excitation energy of the second 
minimum is ~ 1.5 MeV in both models, whereas the superdeformed minimum calculated 
with the SIII interaction is more than 2 MeV higher than in the present calculation. Un- 
less stated otherwise, all calculations presented in this section have been performed in the 
deformed oscillator basis with N s h = 16 oscillator shells. 

The first step in the solution of the modified Hill- Wheeler equation is the construction 
and diagonalization of the norm overlap kernel 

tf(q i ,q j ) = (<f>(q i Mq j )) , (68) 

see Eq. (|60p. Since the GCM basis states are not linearly independent, many of the norm 
overlap kernel eigenvalues are close to zero. This is illustrated in the left panel of Fig. |21 
where we display the eigenvalues nk for four different values of the mesh spacing, ranging 
from Aq = 2b to Aq = 8b. For Aq = 2b the overlaps between neighboring states are 
typically ~ 0.8, and in the corresponding set of 61 eigenvalues we find 13 values smaller 
then 10~ 3 . If the mesh spacing is increased to Aq = 4b and Aq = 6b, the overlaps between 
neighboring states are reduced approximately by factors of two and eight, respectively. A 
further increase of the mesh spacing results in very small overlaps between neighboring states 
(~ 0.05), i.e. basis states become almost orthogonal. Except in a few test cases that will 
be specified explicitly, all calculations in this section have been carried out with the mesh 
spacing Aq = 2b. 

In the next step the GCM basis space is truncated by eliminating those eigenvectors of the 
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norm overlap kernel, which correspond to eigenvalues smaller than a given positive constant 
e n . This is necessary in order to eliminate numerical instabilities in the diagonalization of 
the collective Hamiltonian Eq. (J61|) . In the right panel of Fig. |2l the energies of twelve 
lowest GCM states are plotted as functions of the number of basis states. We notice that 
the spectrum is stable for a broad range of basis dimensions, between 25 (e n = 0.2) and 55 
(e n = 5 • 10 -4 ) vectors. These results can be directly compared with Fig. 4 of Ref. [25j. In 
the following calculations, eigenvectors of the norm overlap kernel with eigenvalues smaller 
than e n = 5 • 10~ 4 are eliminated from the basis. 

In Fig. H2 we plot the energies of fifteen lowest GCM states as functions of the average 
quadrupole moment 



calculated in oscillator bases with 10, 12, 14 and 16 oscillator shells, together with the cor- 
responding mean-field binding energy curves. The GCM ground states are normalized to 
zero energy. In all four cases the average deformation of the ground state is close to the 
minimum of the binding energy curve, and the gain in correlation energy which results from 
configuration mixing is ~ 0.8 MeV. The energies of the ground state and the two first excited 
states basically converge already for a basis with 12 shells. Higher excited states, however, 
contain sizeable admixtures of basis states with larger deformations, and the corresponding 
energy spectrum is sensitive to the number of oscillator shells. 

The GCM states can be analyzed in more detail if one plots their collective wave functions 
gk{q) as functions of the quadrupole moment. In Fig. 0] we display the collective wave 
functions for the first fourteen GCM states in 194 Hg. The vertical dashed line denotes the 
position of the barrier between the main potential well and the superdeformed well. Except 
for the fifth state, the wave functions of the lowest nine states are concentrated in the main 
potential well. The fifth state obviously belongs to the superdeformed minimum, hence its 
energy displays a strong dependence on the number of oscillator shells (see also Fig. |3J). For 
states with k > 10 the wave functions are generally spread over a wide region of deformations, 
both in the main and in the superdeformed well. 

In Fig. El we plot the GCM energy spectra, calculated with 16 oscillator shells, for four 
values of the mesh spacing, ranging from Aq = 2b to Aq = 8b. The corresponding mean- 
field binding energy curves are also included in the figure, and their minima are placed at 




(69) 
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zero energy. Comparing with our standard value of Aq = 2b, we notice that the low-energy 
part of the spectrum is accurately calculated also for Aq = 4b. Increasing the mesh spacing 
to Aq = 6b, accurate energies are obtained only for the two lowest states. With a further 
increase of Ag, the overlaps between neighboring states become so small that there is hardly 
any configuration mixing. The resulting GCM energies are very close to the energies of the 
basis states \<f>(q))- 

Several additional tests, carried out in comparison with the results of Ref. Q|, have 
shown that our implementation of the GCM is numerically stable and, therefore, it can also 
be used for configuration mixing calculations with angular momentum projected states. 



B. Test of angular momentum projection: J Mg 

For a quantitative description of structure phenomena, especially in transitional deformed 
nuclei characterized by a coexistence of spherical and intruder configurations, calculations 
must explicitly include correlations related to restoration of broken symmetries. In partic- 
ular, the rotational energy correction, i.e. the gain in energy obtained by projection on 
states with good angular momentum, can be of the order of 2-4 MeV for the ground state. 
Here we perform several tests of the angular momentum projection for the isotope 32 Mg. 
This nucleus belongs to the island of inversion at N = 20, which is characterized by the 
melting of the neutron shell closure and the predominance of intruder state configurations 
in ground states of neutron-rich systems. The structure of 32 Mg has been the subject of 
numerous experimental and theoretical studies. Several modern theoretical approaches have 
recently been employed in extensive studies of the erosion of the spherical N = 20 shell 
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371 ]. the quantum Monte Carlo shell 



closure in this neutron-rich nucleus: the shell model 

model [3], the angular momentum proj ected generator coordinate method based on the 
non-relativistic Gogny interaction jj, Isflj ]. Although virtually all self-consistent mean-field 
models, non-relativistic as well as relativistic, predict a spherical ground state for 32 Mg, the 
GCM calculation with the Gogny force has shown that the ground state becomes deformed 



as a result of the inclusion of rotational energy correction. Both the excitation ene rgie s 
E(2+) Q and £(4+) QQ, 

as well as the transition probability B(E2, t -> 2+) Q,Q, 
have been measured for 32 Mg. When compared to data from neighboring nuclei, the rela- 
tively low excitation energy of the first excited state E(2^) = 885 keV, the large transition 
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probability B(E2,0^ -> 2^), and the ratio E(4f)/E(2f) = 2.6, indicate that the ground 
state of 32 Mg is deformed. 

In Fig. El we display the mean-field binding energy curves for 32 Mg as functions of the 
quadrupole moment, calculated with the PC-F1 relativistic point-coupling effective inter- 
action. The constrained mean-field equation has been solved self-consistently on a regular 
mesh ranging from q = —2.2b to q = 4.0b, with the mesh spacing Aq = 0.2b. The three 
curves correspond to calculations with N s h = 8, 10 and 12 major oscillator shells. For such 
a light system and for this range of deformations, it appears that already 10 oscillator shells 
are sufficient to obtain a reasonably converged mean-field binding energy curve. In the 
following calculations we expand the Dirac spinors in the axially deformed oscillator basis 
with N s h = 10 major shells. This choice is also supported by the results of Ref. jj], where 
correlations beyond the mean-field approximation have been studied in the framework of the 
angular momentum projected GCM with the Gogny force. In addition to a spherical ground- 
state, the PC-F1 binding energy curves display a prolate deformed shoulder at q — 1.5b, 
at an excitation energy of ~ 3.5 MeV above the ground state. The binding energy curve 
calculated with the Gogny force is similar (see Fig. 6 of Ref. [i]), but the prolate shoulder 
is somewhat more pronounced, and is located only « 1.9 MeV above the spherical ground 
state. Of course, if the shoulder is too high above the spherical ground state, correlations 
related to the restoration of rotational symmetry and quadrupole fluctuations might not 
be strong enough to deform the nucleus. The different predictions for the location of the 
shoulder can be related to the single-particle levels calculated with the PC-F1 interaction, 
displayed in Fig. and with the Gogny interaction (Fig. 5 of Ref. |4|). In these figures 
the eigenvalues of the corresponding mean-field Hamiltonians are plotted as functions of the 
quadrupole deformation. The ratio between the neutron spherical gap (7.2 MeV), and the 
gap at deformation q = 1.5b (2.9 MeV), is ~ 2.5 for the PC-F1 interaction, whereas the 
Gogny force gives a much smaller value for this ratio « 1.8. This leads to a more pronounced 
prolate shoulder at lower excitation energy. 

The essential step in the procedure of angular momentum projection is the evaluation of 
the projected norm overlap kernel 

Af J (q,q) = (<j>(q)\P&\<j>(q)) = (2 J + 1) 1 + t*¥- j*'* sin /3d&(J3)n(q; , (70) 
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where 

n(g;^) = (0(g)|e-^ |0(g)) . (71) 

In several studies 135. 39 



it has been shown that the ansatz 

n app (q; (3) = e-^ J >^ , (72) 

presents an excellent approximation for the function n(q; f3), both at small and large defor- 
mations. The expectation value ( Jy ) as a function of the quadrupole moment is plotted in 
Fig. |HJ This curve is in agreement with the one obtained with the Gogny interaction (see the 
right panel in Fig. 5 of Ref. I n Fig. El we display the function n(q; (3) for several val- 

ues of the quadrupole moment. The solid curves correspond to the approximate expression 
Eq. (j72j) . whereas dots denote values obtained with the exact calculation. The comparison 
between the exact and approximate results provides a very useful test of the numerical pro- 
cedure used in angular momentum projection. The projected norm overlap kernels, shown 



in Fig. El for the four lowest angular momenta, can 



De compared with those obtained using 



4] 



i. We notice that the spherical con- 



the Gogny effective interaction (see Fig. 7 of Ref. 
figuration is a pure + state (A / " J=0 (0, 0) = 1). The maxima of the projected norm overlap 
kernels for higher angular momenta are correspondingly shifted to larger deformations. 

In Fig. ^2 the energies of the angular momentum projected states are analyzed. At this 
stage we do not consider configuration mixing yet, and the projected energy of the \<f>(q)) 
state reads 

The angular momentum projected energy curves for J 71 " = + ,2 + ,4 + ,6 + , and 8 + are plot- 
ted, together with the corresponding mean-field binding energy curves, as functions of the 
quadrupole deformation. The curves obtained from solutions in axially deformed oscillator 
bases with N s h = 8, 10 and 12 major shells are almost identical. Since the spherical configu- 
ration is already a pure + state, there is no energy gain for J 71 " = + at q = Ob. Notice that 
the spherical point q = Ob is not included in plots of E J (q) for J > 2. Namely, for J ^ the 
quantities H J (0, 0) and A/" J (0, 0) are so small, that their ratio Eq. ()73|) cannot be determined 
accurately. For higher values of the angular momentum (J 71 = 6 + , 8 + in Fig. ITT|l several 
additional configurations close to the spherical point are also characterized by very small 
values of the projected norm overlap kernel. These configurations can be safely omitted 
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from the projected energy curves, because on the one hand the angular momentum projec- 
tion becomes inaccurate at these points, and on the other hand the corresponding angular 
momentum projected states would not play any role in configuration mixing calculations. 

Is is interesting to compare the projected energy curves with those obtained using the 
Gogny effective interaction (see Fig. 6 of Ref. |4J). The principal difference is seen already 
for the J = + projected energy. The PC-F1 interaction predicts two almost degenerate 
minima at small oblate and prolate deformations. The occurrence of degenerate oblate and 
prolate minima, symmetrical with respect to the spherical configuration, is a feature com- 
mon to all nuclei for which the mean-field calculation predicts a spherical ground state [39| . 
As compared to the mean-field energy, the prolate deformed shoulder is more pronounced for 
the E J=0 (q) curve, and its excitation energy has been lowered from 3.5 MeV to 1.2 MeV by 
angular momentum projection. On the other hand, at the mean-field level the Gogny inter- 
action predicts a more pronounced shoulder, only « 1.9 MeV above the spherical minimum. 
With angular momentum projection the shoulder becomes the absolute minimum of the 
J = projected energy curve. Therefore, the inclusion of the rotational energy correction 
leads to a deformed ground state in 32 Mg, when calculated with the Gogny interaction. In 
addition, the degenerate oblate and prolate minima, symmetrical with respect to the q = 0, 
are predicted at slightly higher excitation energy. In the present calculation with the PC-F1 
interaction, the gain in rotational energy is too small to deform the ground state of 32 Mg. 
The rotational energy correction -Erec ; i-e. the difference between the mean-field and the 
.P = + projected energy curves is plotted in Fig. H21 -Erec is zero for the spherical intrin- 
sic state, and generally it increases rather steeply for small deformations (see also Fig. ITTf . 
Our result for -Erec is very similar to the curve obtained from the Gogny mean-field poten- 
tial energy (see Fig. 6 of Ref. j^]])- This means that the deviation between the J 77 = + 
projected energy curves can indeed be attributed to the difference between the PC-F1 and 
Gogny interactions on the mean-field level. 

C. Angular momentum projection and configuration mixing: 32 Mg 

As a final test of our implementation of the GCM for relativistic mean-field models, 
we have performed configuration mixing calculations of the angular momentum projected 
intrinsic states for 32 Mg. The solution of the Hill- Wheeler matrix equation (|fi2*|) . with the 
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collective Hamiltonian Eq. (|6ip. determines both the ground-state energy, and the energies 
of excited states, for each value of the angular momentum J. The collective wave functions 
5a(<?)> an d the weight functions are calculated from the norm overlap eigenfunctions 

Eqs. (|63j) and (J64j) . respectively. As we have shown in Fig. fo r J > 2 several points 
on the energy surfaces close to the q = correspond to configurations with very small 
values of the projected norm kernel. Since the numerical evaluation of the norm overlap and 
Hamiltonian kernels is not accurate in such cases, we have excluded from the configuration 
mixing calculation all those intrinsic configurations for which Af J (q, q) < 0.001 . 

The energies and the average quadrupole moments Eq. (jfi9j) of the two lowest GCM 
states for each angular momentum are displayed in Fig.^3 together with the corresponding 
projected energy curves. The spectrum can be compared with the available data, and 
with the angular momentum projected GCM results obtained using the Gogny effective 
interaction (see Fig. 7 of Ref. |39|). Configuration mixing between the two essentially 
degenerate oblate and prolate minima of the J = + energy curve, symmetrical with respect 
to q — 0, results in the almost spherical ground state 0^. When calculated with the Gogny 
interaction, on the other hand, the ground state is prolate deformed. The relatively large 
deformation is a result of a fine balance between the energy correction associated with the 
restoration of rotational symmetry (favors larger deformation), and the correlations induced 
by quadrupole fluctuations (mixing between oblate and prolate configurations reduces the 
deformation of the lowest + state). 

The excitation energies of the 2f, 4f and 6+ GCM states are included in TablelH together 
with the corresponding energies obtained with the Gogny force, and the available experimen- 
tal excitation energies. Obviously, the PC-F1 interaction predicts yrast states at excitation 
energies that are too high, compared with the Gogny interaction, or with the experimental 
values. In Table ITT1 we display the spectroscopic quadrupole moments of the 2f, 4+ and 
6^ GCM states, for the PC-F1 and Gogny effective interactions. Although comparable 
in size, the quadrupole moments calculated with the PC-F1 interaction are systematically 
smaller. This is, of course, consistent with the lower excitation energies predicted by the 
Gogny force. Since the ground state is almost spherical, the calculated transition probability 
B(E2; 0± — > 2f) = 15.5 e 2 fm A , is far too small when compared to the experimental value 
447(57) e 2 fm 4 ) 4^, or to the value obtained with the Gogny interaction (395 e 2 /m 4 ) 



39]. 
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The differences in the spectra predicted by the PC-F1 and Gogny interactions originate in 
the deviation of the corresponding mean-field binding energy curves or, more precisely, in the 
different neutron single-particle levels (Nilsson diagrams) calculated with the two effective 
interactions. Because of the large spherical gap predicted by the PC-F1 interaction, the 
magic number N = 20 persists even in such a neutron-rich system, and 32 Mg exhibits 
structure properties typical for other magic nuclei, e.g. 48 Ca [39]. This is further illustrated 
in Fig. HU where we display the amplitudes of the collective wave functions \gi(q)\ 2 for 
the two lowest GCM states of each angular momentum, together with the corresponding 
projected energy curves. For instance, \g\{q)\ 2 obviously reflects a configuration mixing 
of the prolate and oblate minima with almost equal weights, resulting in a ground state 
with an average quadrupole moment close to zero. The rotational energy correction for 
the ground state, i.e. the energy gain from angular momentum projection, is ~ 1 MeV. 
Configuration mixing provides an additional gain of 0.3 MeV. Both values are in agreement 
with the corresponding quantities calculated for the magic 48 Ca nucleus The amplitudes 
\g{ (q)\ 2 for J = 2, 4 and 6, are localized in the prolate wells of the corresponding projected 
energy curves and, therefore, the average quadrupole moments of the states 2+, 4^ and 6^ 
are close to the prolate minima. The collective wave functions g£(q), for J = 2 ,4 and 6, 
correspond to a band based on the (3 vibrational state 0^. 

IV. SUMMARY AND OUTLOOK 

The framework of self-consistent relativistic mean-field (RMF) models has been very 
successfully employed in analyses of a variety of nuclear structure phenomena, not only in 
nuclei along the valley of /3-stability, but also in exotic nuclei with extreme isospin values 
and close to the particle drip lines. Applications have reached a level of sophistication 
and accuracy comparable to the non-relativistic Hartree-Fock (Bogoliubov) approach based 
on Skyrme or Gogny effective interactions. Although mean-field and pairing correlations 
are treated very carefully in modern RMF models, additional correlations, related to the 
restoration of broken symmetries and to fluctuations, have either been neglected or taken into 
account in an implicit way. In this work we have introduced a model in which restoration of 
rotational symmetry and fluctuations of the quadrupole deformation are explicitly included 
in the relativistic framework. 
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In the specific model which has been developed in this work, the generator coordinate 
method (GCM) is employed to perform configuration mixing calculations of angular mo- 
mentum projected wave functions, calculated in a relativistic point-coupling model. The 
geometry has been restricted to axially symmetric shapes, and the mass quadrupole mo- 
ment is used as the generating coordinate. The intrinsic wave functions are generated from 
the solutions of the constrained relativistic mean-field + BCS equations in an axially de- 
formed oscillator basis. 

In order to test our implementation of the GCM and angular momentum projection, a 
number of illustrative calculations have been carried out for the nuclei 194 Hg and 32 Mg. The 
PC-F1 parameter set ^6| has been used for the effective point-coupling Lagrangian, and the 
effective interaction in the particle-particle channel has been approximated by a density- 
independent 5-force. The test of the generator coordinate method has been performed in a 
study of quadrupole dynamics in the nucleus 194 Hg, and the results have been compared with 
the classical analysis of the GCM in the investigation of shape isomerism in 194 Hg by Bonche 
et al. (2^1, based on the non-relativistic constrained Hartree-Fock plus BCS model with the 
SIII Skyrme effective interaction. Angular momentum projection and, finally, configuration 
mixing of angular momentum projected states, have been tested in the example of the 
neutron-rich nucleus 32 Mg, in comparison with results obtained with the angular momentum 
projected generator coordinate method based on the non-relativistic HFB with the Gogny 
interaction jj, The tests have been very successful, and the results obtained for the 
binding energy curves, projected energy curves, rotational energy corrections, ground and 
low-lying excited states, and collective wave functions for 194 Hg and 32 Mg, are generally in 
very good agreement with the predictions of GCM calculations based on non-relativistic 
Skyrme and Gogny interactions, respectively. 

The choice of the PC-F1 relativistic effective interaction, however, does not lead to a de- 
formed solution for the ground state of 32 Mg, even after the inclusion of the rotational energy 
correction. This result is in contrast with available data, and with the configuration mixing 
calculation of angular momentum projected configurations based on the Gogny interaction. 
The different predictions for the ground state of 32 Mg can be related to the corresponding 
mean-field binding energy curves and, more specifically, to the different results for the size of 
the spherical N = 20 neutron gap, obtained with the Gogny and PC-F1 interactions. Even 
though the spherical ground state of 32 Mg, predicted by the PC-F1 effective interaction, is 
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not crucial in the context of the present analysis, it points to an important problem. Namely, 
the choice of effective interactions to be used in self-consistent calculations that go beyond 
the mean-field approximation and explictly include correlations, such as those considered in 
the present work. Virtually all global effective interactions have been adjusted to data, e.g. 
masses and radii, which already include correlations. On the other hand, those correlations 
that we wish to treat explicitly, should not be included in the effective interaction in an 
implicit way. The solution is to adjust global effective interactions to pseudodata, obtained 
by subtracting correlation effects from experimental masses and radii. Approximate meth- 
ods for the calculation of correlations have recently been developed [35] , that will enable a 
systematic evaluation of correlation energies for the nuclear mass table. 

Before proceeding with realistic applications of the model introduced in this work, our first 
task is to adjust a new global effective point-coupling interaction which will not implicitly 
contain rotational energy corrections and quadrupole fluctuation correlations. Further de- 
velopments will include the treatment of pairing fluctuations by particle number projection, 
the use of different generating coordinates for the neutron and proton density distributions, 
the description of non-axial shapes, and the extension to odd nuclei. 
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FIG. 1: The binding energy curves for Hg, as functions of the mass quadrupole moment, calcu- 
lated by expanding the Dirac spinors in 10, 12, 14, 16 and 18 oscillator shells. 

FIG. 2: The eigenvalues of the norm overlap kernel, calculated using four different values for the 
mesh spacing: Aq = 2,4,6 and 8 b (left panel). The energies of the twelve lowest GCM states 
plotted as functions of the dimension of the GCM basis (right panel). 

FIG. 4: GCM collective wave functions gk(q) Eq. (|57jl for the lowest fourteen states in 194 Hg. The 
vertical dashed line denotes the position of the barrier separating the main and the superdeformed 
potential wells. 

FIG. 5: The mean-field binding energy curves for 194 Hg, together with the energies and average 
quadrupole moments of fifteen lowest GCM states. Calculations have been performed using four 
values of the mesh spacing: Aq = 2, 4, 6 and 8 b. Zero energy is placed at the position of the 
minimum of the binding energy curve. 

FIG. 6: The binding energy curves for 32 Mg, calculated from the constrained solutions of the self- 
consistent relativistic mean-field equations in axially deformed oscillator bases with 8, 10 and 12 
major shells. 

FIG. 7: The neutron (left panel) and proton (right panel) single-particle levels for 32 Mg, as functions 
of the mass quadrupole moment. The thick dashed curve denotes the position of the Fermi energy. 

FIG. 8: The expectation value (J 2 ) for 32 Mg, as a function of the mass quadrupole moment. 

FIG. 9: A comparison between the exact values of the function n(q; [3) Eq. (|71|) (dots) and the 
approximate expression Eq. (|72[) (curves), for several values of the mass quadrupole moment. 

FIG. 3: The energies of fifteen lowest GCM states in 194 Hg, plotted as functions of the average 
quadrupole moment, together with the corresponding mean-field binding energy curves. The four 
panels correspond to calculations in oscillator bases with 10, 12, 14 and 16 major oscillator shells. 
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FIG. 10: Projected norm overlap kernel Af J (q, q) as a function of the mass quadrupole moment for 
32 Mg. 



FIG. 11: Angular momentum projected (J w = + ,2 + ,4 + ,6 + and 8 + ) potential energy curves for 
32 Mg, as functions of the mass quadrupole moment. The mean-field energies are also included 
(thick dotted curves). The three panels correspond to solutions in axially deformed oscillator bases 
with 8, 10 and 12 major shells. 

FIG. 12: Rotational energy correction as a function of the mass quadrupole moment for 32 Mg. 

FIG. 13: The energies and the average quadrupole moments of the two lowest GCM states for each 
angular momentum in 32 Mg, together with the corresponding projected energy curves. 

FIG. 14: Squares of the collective wave functions \gi{q)\ 2 of the two lowest GCM states for each 
value of the angular momentum in 32 Mg, together with the corresponding projected energy curves. 

TABLE I: The excitation energies (in MeV) of the 2f, and 6^ GCM states for 32 Mg. 





E (PC-F1) 


E (Gogny) E (exp.) 




2.04 


1.4 0.885 




4.42 


3.6 1.437 




7.41 


5.5 



TABLE II: The spectroscopic quadrupole moments (in e fm 2 ) of the 2f, if and 6^ GCM states 
for 32 Mg. 





Q spec (p C _ F1 ) 


QT C (Gogny) 




-17.51 


-19.15 




-19.28 


-26.31 




-21.19 


-30.09 
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